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Abstract 

The improved quark mass density-dependent (IQMDD) model, which has been successfully used 
to describe the properties of both infinite nuclear matter and finite nuclei, is applied to investigate 
the properties of quark deconfinement phase transition. By using the finite-temperature quantum 
field theory, we calculate the finite temperature effective potential and extend the IQMDD model 
to finite temperature and finite nuclear matter density. The critical temperature and the critical 
density of nuclear matter are given and the QCD phase diagram is addressed. It is shown that 
this model can not only describe the saturation properties of nuclear matter, but also explain the 
quark deconfinement phase transition successfully. 
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Owing to the nonperturbative nature of quantum chromodynamics (QCD) in low energy 
regions, it is very difficult to study nuclear system by using QCD directly. Phenomenological 
models reflecting the characteristic of the strong interaction are widely used in the studying 
of the properties of hadrons and nuclear matter. The quark-meson coupling (QMC) model 
suggested by Guichon [l| is a famous hybrid quark meson model, which can describe the 
saturation properties of nuclear matter and many other properties of finite nuclei successfully 
[2,3]. In this model, the nuclear system was suggested as a collection of MIT bag, vector 
w-meson and scalar cx-meson, and the interactions between quarks and mesons are limited 
within the MIT bag regions due to the quark cannot escape from the MIT bag. 

Although the QMC model is successful for describing the physical properties of a nuclear 
system, two shortcomings arise when one uses this model to discuss some physical topics. 
First, the MIT bag model is a permanent quark confinement model and its boundary condi- 
tion cannot be destroyed by temperature and nuclear matter density. The Second difficulty 
arises from the MIT boundary condition. If we hope to do the nuclear many-body calcula- 
tions beyond mean field approximation (MFA) by quantum field theory, it is essential to find 
the free propagators of quark, a meson and u meson, respectively. However, the constraint 
of the MIT bag boundary condition presents obstacles to get the corresponding propagators 
in free space. Because the interactions between quarks and mesons are limited within the 
bag regions, and multireflection of quarks and mesons by the boundary must be taken into 



account {4]. 

In order to keep the quark confined property and give up the MIT boundary condition, a 
new quark-meson coupling model based on a quark mass density-dependent (QMDD) model 
was presented in our previous papers [5,6]. According to the QMDD model, the masses of 
u, d and s quarks (and the corresponding antiquarks) satisfy: 

m q = - — (q = u,d,u,d) } (1) 
6n B 

B 

m St s = m s0 + - — , (2) 
3n B 

where Ub = \{ n u + n d + n s ) is the baryon number density, with n u ,n d ,n s representing 
the density of the u quark, d quark, and s quark, m s o is the current mass of the strange 
quark, and B is the vacuum energy density inside the bag. As was explained in Refs. [7|, 
the ansatz Eqs. (1) and (2) corresponds to a quark confinement hypothesis because when 



V — > oo, ub — > and m q — > oo, it prevents the quark from going to infinity or to very 
large regions. The large volume means that the baryon density is small. This mechanism of 
confinement can be mimicked through the requirement that the mass of an isolated quark 
become infinitely large so that the vacuum is unable to support it [8j. This is just the 
physical picture given by Eqs. (1) and (2). In fact, this confinement mechanism is very 
similar to that of the MIT bag model But the advantage of the QMDD model is that it 
does not need to introduce a quark confined boundary condition like that of the MIT bag 
model. 

But the QMDD model still has two shortcomings: (1) it is still an ideal quark gas model. 
No interactions between quarks exist except for a confinement ansatz, equation (1) and (2). 
(2) it cannot describe the quark deconfinement phase transition and predict a correct phase 
diagram as that given by lattice QCD. The reason is that the temperature T tends to infinite 
when density n B — > 0. This result can easily be understood if we notice the basic hypothesis 
Eqs. (1) and (2) of the QMDD model, the quark masses are divergent when % — > 0. To 
excite an infinite weight particle, one must prepare to pay the price for infinite energy, i.e., 
infinite temperature. It means that the confinement in the QMDD model is still permanent. 
To overcome the difficulty (2), we have introduced a new ansatz that the vacuum density B 
is a function of temperature T [101 ]: 

B(T) = B [l - (T/T c ) 2 ],0 < T < T c , (3) 
B(T) = 0,T>T c . (4) 

and extended the QMDD model to a quark mass density- and temperature- dependent 



(QMDTD) model 10]. This QMDTD model has been employed to discuss the properties of 



strange quark matter [ll| . the dibaryon system [12[] successfully. 

To improve the shortcoming (1) of the QMDD model and to avoid the ad hoc ansatz of 
QMDTD model, we introduced the w-meson and u-meson in the QMDD model to mimic the 
repulsive and attractive interactions between quarks in Refs. [5,6]. Following the treatment 
of Friedberg and Lee [13,14], the interaction between quarks and the nonlinear a field forms 
a Friedberg-Lee bag [§]. In Ref. [5|, we found the wave functions of the quark ground state 
and the lowest one-quark excited states. By using these wave functions, we calculated many 
physical quantities such as root-mean-square radius, the magnetic moment of nucleon to 
compare with experiments and come to a conclusion that this IQMDD model is successful 



to explain the properties of nucleoli. Then under mean field approximation, we employ 
the IQMDD to investigate the physical properties of nuclear matter. Instead of the MIT 
bag, after introducing the nonlinear interaction of cr-mesons and qqu coupling, we construct 
a Friedberg- Lee soliton bag in nuclear system. The quark and a- meson coupling and 
the quark and uj— meson coupling are introduced to mimic the attractive and the repulsive 
interactions between quarks in this model. It has been shown that our model can successfully 
describe the saturation properties, the equation of state, the compressibility and the effective 
nucleon mass of symmetric nuclear matter [6j and finite nuclei 15]. 

We hope to emphasize that there are two basic differences between the IQMDD model 
and the usual QMC model. First, unlike in the QMC model, we do not need an MIT bag 
for the nucleon in the IQMDD model. The constraint of the MIT bag boundary condition 
disappears in our formulas because the quark confinement mechanism has been established 
in Eqs. (1) and (2). Second, the interaction between the quark and the scalar meson is 
limited in the bag regions for the QMC model, while for the IQMDD model this interaction 
is extended to the whole free space. So it is convenient to write down the free propagators 
of quarks and mesons. 

In Ref. [16,17], we extended this model to finite temperature and studied its soliton 
solution by means of the finite temperature quantum field theory. The critical temperature 
of quark deconfinement Tq and the function of temperature-dependent bag constant B(T) 
are found as an output. But the the deconfinement properties of the IQMDD model at 
finite density of nuclear matter and the QCD phase diagram have not yet been studied. We 
would like to study these problems in this paper. We will show that we can get reasonable 
quark deconfinement critical temperature and critical density of nuclear matter. Moreover 
the QCD phase diagram can also be given in the IQMDD model phenomenologically. 

Before giving the main formulas of the IQMDD model, we should emphasize that the 
quark deconfinement phase transition and QCD phase diagram have attracted a lot of in- 
terest for many years, many theoretical efforts have been recently and less recently carried 
out, on the study of the hadronic phase structure and the deconfinement phase transition 
by using the Nambu-Jona-Lasinio model and some its extensions in the framework of finite 
temperature field theory [19-30]. Among them, some perform a very detailed investigation 
of the QCD phase structure and other ones give precise numerical estimates for the decon- 
finement temperature and size of hadrons. These publications is very interesting and also 



4 



give us lots of important information for the QCD phase structure. 
The effective Lagrangian density of the IQMDD model is given by 



C = #7^ 



m„ 



9> 



(5) 



where Fn U = dn0o u — d u Uu, ip represents the quark field, m q 



B 

3n B 



is the mass of u(d) quark, 



the a and u fields are not dependent on time, is the coupling constant between the quark 
field ip and the scalar meson field a, is the coupling constant between the quark field ip 
and the vector meson field uj^ U(a) is the self interaction potential for a field. The potential 
field U(a) is chosen as [14 1 

b 



b > 3ac. 



The condition ([7]) ensures that the absolute minimum of U(a) is at a 
potential U(a) has two minima: one is the absolute minimum a v , 



2c 



1 + 



1 - 



8ac 
362 



(6) 
(7) 

<7 V ^ 0. The 

(8) 



it corresponds to the physical vacuum, and the other is at <j = 0, it represents a metastable 
local false vacuum. We take U(a v ) = and the bag constant B can be expressed as 



d o bo C A 



(9) 



In order to study the deconfinement phase transition, we turn to extend IQMDD model to 
finite temperature. The appropriate framework is finite temperature quantum field theory. 
The finite temperature effective potential plays a central role within this framework. Under 
the mean field approximation, the oj meson field operator can be replaced by its expectation 
value in symmetric nuclear matter 18], 



A 3g " 

6 ^~^tPb- 



(10) 



Using the method of Dolan and Jackiw [31[, up to one- loop approximation, the effective 
potential reads: 



V{a- (3- VJ) = U(a) + V B (a; (3) + V F (a; V u ), 



(11) 



where 

V u = gftD = 3 9 -^ P b (12) 
mj 

is the contribution of u- field, (3 is the inverse of the temperature, ps is the density of nuclear 
matter and 



2tt 



l — J™ dxx 2 \n (l - e-V^) , (13) 

V F (a; 0; VJ) = -12 £ ^ cte 2 ln (l + e -<V^^-/%>) , (14) 

where the minus sign of Eq. (14) is the consequence of Fermi-Dirac statistics. The degenerate 
factor 12 comes from: 2(particle and antiparticle), 2 (spin), 3(color). m a and m q are the 
effective masses of the scalar field a and the quark field respectively which can be found in 
Ref. [16,17]. 

We see from Eqs. (11-14) that the scalarlike interaction g^J + aJ gives contribution for 
effective masses of quark and a meson and then forms a confined soliton bag, while the 
vectorlike interaction g^ip + 'y fl u fl il ! gives contribution to an effective chemical potential of 
quarks. In fact, this finite temperature effective potential for FL model had been calculated 
by many others authors [16,17]. 

In the FL bag model, the finite-temperature vacuum energy density B(T) is defined as 

BJ3; VJ = V(a ; (3; VJ) - V{a v ; J VJ. (15) 

It is the difference between the values in the perturbative false vacuum state and the values 
in the physical real vacuum state of the finite temperature effective potential. At critical 
density of nuclear matter and critical temperature of quark deconfinement phase transition, 
B equal to zero: B(T,p B ) = 0. 

Before numerical calculation, let us discuss the parameters in IQMDD model. First, we 
choose rrij = 783 MeV, m p = 770 MeV and m a = 509 MeV as that of Ref. fl. Fixing the 
nucleon mass = 939 MeV, we get B = 174 MeV fm -3 . Obviously, the behaviors at the 
saturation point must be explained for a successful model. It reveals that nuclear matter 
saturates at a density po = 0.15 fm -3 with a binding energy per particle E/A— —15 MeV at 
zero temperature, and the compression constant to be about K(p Q ) = 210 MeV. Therefore 
we fixed g* = 2.44, gfc = 4.67, b = -1460 MeV to explain above data 



Now we turn to investigate the effective potential V(a; (3; Vj). The effective potential at 
finite temperature and finite nuclear matter density can be obtained by numerical calcula- 
tions using the set of Eqs. (11-14). Then it is convenient to investigate the temperature 
and the nuclear matter density dependence of the bag constant. The bag constant in the 
IQMDD model is defined as Eqs. (15). The upper panel of Fig. 1 shows the temperature 
dependent behavior of the bag constant at various nuclear matter densities. One can see 
clearly that the bag constant decreases continuously with increasing temperature at fixed 
nuclear matter density. At critical temperature Tc, the bag constant equals to zero and 
the deconfinement phase transition begins to occur. The lower panel of Fig. 1 shows the 
nuclear matter density dependence of the bag constant at different temperatures. Similar 
to the upper panel of Fig. 1, one can see the deconfinement phase transition begin to take 
place at a critical density for a given value of temperature. 

The critical behavior of nuclear matter in the T - p B phase diagram calculated using the 
IQMDD model is described in Fig. 2. The critical line corresponds to the phase transition 
between confined and deconfined phases. In this nuclear model, our calculation shows that 
the critical density of nuclear matter for quark deconfinement is about ten times the sat- 
uration density of the nuclear matter. The phase diagram shown in Fig. 2, based on the 
phenomenological Friedberg-Lee model, predicts a first-order deconfinement phase transi- 
tion for the full phase diagram. This is due to at the critical T and p B the two vacuum 
states, namely the perturbative vacuum state and the physical vacuum state, appear in the 
same energy. This result differs from the predictions based on lattice gauge theory where 
more complicated phase diagram on the T and p B plane has been obtained, giving more 
fruitful behaviors of the QCD phase [20,21]. For example, at finite temperature and zero 
chemical potential, the phase transition could be of a second order; moreover, it is predicted 
that there may exist a critical point along the critical line. 

In the end of paper, we have provided a comparison of the results obtained from the 
IQMDD model with the recent results obtained in the framework of effective four fermion 
interaction theories. The recent efforts [21-30] have made very fruitful progress on the QCD 
phase structure including finite-size hadron and magnetic effects. Their theoretical frame 
of field theory is the imaginary-time temperature field theory, same as the IQMDD model. 
As the IQMDD model has lack of chirality, we can not predict critical point along the 
critical line. However, we should emphasize that the IQMDD model can not only explain 
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the quark deconfiriement phase transition, but also describe the saturation properties of the 
nuclear matter successfully. In the IQMDD model, the symmetry nuclear matter is treated 
as uniform and unlimited, the vectorlike interaction g^ip^^uj^ip between quarks and u> 
meson gives contribution to an effective chemical potential of quarks, and the expected 
value of u meson is connected with the baryon density. Our numerical results shows that 
the deconfming baryon density at zero temperature is about 1.63 fm -3 which is very close 
to the value predicted by the lattice QCD theory. Meanwhile, the deconfming temperature 
at zero baryon density is 139 MeV which is also reasonable compared to 170 MeV given 
by lattice QCD theory |33|. Recently, the Science journal |34j report that the deconfming 
critical temperature for the QCD phase transition at zero baryon density is about 175 MeV. 

In conclusion, we investigate the quark deconfmement phase transition at finite temper- 
ature and finite density of symmetric nuclear matter in the IQMDD model. We would like 
to emphasize that this is the basic advantage of IQMDD model, because the saturation 
properties can not only be explained by IQMDD model but also by QMC model. The 
reason for quark deconfmement can be explained by IQMDD model is that MIT boundary 
constraint has been dropped out and interactions between quark and mesons have been 
extended to the whole space. Instead of the MIT bag in QMC model, a FL soliton 
bag is introduced in the IQMDD model, which makes it possible to discuss the quark 
deconfmement phase transition. The spontaneously breaking symmetry of nonlinear a field 
is restored and the soliton bag will disappear at critical condition. 

The author wish to thank Professor Ru-keng Su for useful correspondences. This work 
is supported by the National Natural Science Foundation of China (Grants 11105072 and 
11035009). 
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FIG. 1: (upper panel) Bag constant vs temperature at different nuclear matter density. A: ps = 0.0 
fm -3 , B: pb = 0.8 fm~ 3 , C: pb = 1.2 fm~ 3 , D: ps = 1.5 fm~ 3 ; (lower panel) Bag constant vs nuclear 
matter density at different temperature. A: T = MeV, B: T = 60 MeV, C: T = 120 MeV. 
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FIG. 2: Phase diagram of deconflnement in the IQMDD model. 



12 



